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ABSTRACT 

We investigate the possibility that gamma-ray bursts orig- 
inate from quakes deep in the solid crust of a neutron star. 
Seismic waves are radiated if shear stress is relieved by 
brittle fracture. However they cannot propagate directly 
to the surface but are temporarily trapped below a re- 
flecting layer. The shaking of the stellar surface couples 
the seismic waves to Alfven waves which propagate out 
into the magnetosphere. The crust-magnetosphere trans- 
mission coefficient strongly increases with wave frequency 
and magnetic field strength. Alfven wave luminosities suf- 
ficient to power galactic gamma-ray bursts are possible if 
magnetic fields > 10 n G cover at least part of the stellar 
surface. As the Alfven waves propagate out into the low 
density magnetosphere, they become increasingly charge 
starved, thereby accelerating particles to relativistic ener- 
gies. 
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1. INTRODUCTION 


discussion of this model can be found in Ref. 8. 

2. QUAKE ENERGETICS 

The energy released in a typical gamma-ray burst, as- 
suming isotropic emission, is 

„ 10 38 ( F \ ( D \ 1 2 

Eb *~ VlO-«ergscm-^(lk^J ^ 

where 77 specifies the fraction of energy released which is 
converted into gamma-rays, F is the observed fluence, and 
D is the distance. 

We model the crust as plane parallel, chemically homo- 
geneous, and subject to a constant Newtonian gravita- 
tional acceleration, g = — gz , where 2 is a unit vector 
in the upward direction. In numerical expressions we take 
g = 10 14 cms -2 . Above neutron drip, the pressure is due 
to degenerate electrons. Integrating the equation of hydro- 
static equilibrium, we obtain the density profile 


Fifteen years after their discovery (Ref. 1), gamma- 

ray bursts remain a major unsolved problem of high energy 

astrophysics. The bursts lack counterparts at other wave- 
lengths, so the source objects are still unidentified. Neu- 
tron stars have long been viewed as being the most plau- 
sible sources. Individual bursts usually last for a few sec- 
onds, with rise times as short as several milliseconds (Ref. 
2). The relative lack of low-energy X-rays, L K /L y ~ 0.02, 
is an important clue to understanding the nature of the 
gamma-ray emission region. This X-ray paucity constraint 
(Ref. 3) is unlikely to be satisfied in regions of high density 
or near the neutron star surface, where substantial ther- 
mal reprocessing of gamma-rays into X-rays would occur. 
It suggests emission from a region of low- density plasma 
well out in the magnetosphere. 

The basic idea of starquake models for gamma-ray 
bursts is quite simple (Refs. 4-7). Elastic energy released 
in a crust quake excites oscillations of the magnetic field 
frozen in the surface, and the induced electric field accel- 
erates high-energy particles which in turn radiate gamma- 
rays. This work is devoted to assessing the merits of star- 
quakes as the root cause of gamma-ray bursts. A detailed 
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+ 2gm e 



( 2 ) 


where m u is the atomic mass unit, p e is the mean molecular 
weight per electron, and all other symbols have their usual 
meanings. We adopt an intermediate value of /i e = 2.5 
when presenting numerical results. 

The ions in the solid crust are arranged in a Coulomb 
lattice whose shear modulus is given by 

n = 0.295ZVn 8/ \ (3) 

where is the ion number density and Z is the atomic 
number. An average value of Z — 32 is adopted here. The 
elastic energy released in a quake may be expressed as 


Eq ~ fi( z Q)e 2 yi ' ld d 3 * } (4) 

where e y j is the yield strain at which the crust cracks, 
d 2 is the area of the fault plane, and zq is the depth. 
Equations (3) and (4) together require: 
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d ~ 3 x 10 4 


Eq \ 1/3 

10 33 ergs ) V 10- J / 


cm, 


( 5 ) 


where, for convenience, the quake is assumed to take place 
at neutron drip. 

The characteristic frequency, i/ 0 , of elastic waves emit- 
ted by the quake is simply the speed at which the frac- 
ture propagates divided by d . Since cracks typically prop- 
agate at a significant fraction of the shear wave speed, 
v t = (/x/p) 1 / 2 , we find 



The seismic energy density is partitioned into shear 
and compressional waves in proportion to the inverse sixth 
power of the ratio of their propagation speeds. Most of 
the energy is emitted in shear waves since their speed is 
typically a few times smaller than that of compressional 
waves. 


3. MAGNETOELASTODYNAMICS 
The magnetic energy density dominates the rigidity 
in the outer layers of the crust of magnetic neutron stars 
and exceeds the rest mass energy density in the magneto- 
sphere. To include the effects of the magnetic field, B, on 
wave propagation, we add the Maxwell stress to the equa- 
tions of elasto dynamics. For simplicity, we consider the 
unperturbed magnetic field to be uniform and we treat 
the surface layers as solid because, with our scalings, the 
magnetic field completely dominates the stress in the ocean 
provided the surface temperature is < 3 X 10 6 K. This is 
probably true for old neutron stars. 

Linearizing the equation of motion and the continuity 
equation about a static equilibrium, we obtain 


+ -Sj x B + 6pg - VS Pt (7) 
di 2 c 

6p = - V-(pO, ( 8 ) 

respectively. Here, 6 denotes an Eulerian perturbation, £ 
is the displacement of an element of material from its equi- 
librium position, and Sj is the perturbed current density. 
The pressure due to the degenerate electrons, p, has been 
written separately from the elastic stress associated with 
the deformation of the Coulomb lattice, <x, The compo- 
nents of the perturbed elastic stress tensor are related to 
the components of the gradient of £ by 


= ( k - f) ***** + ** (Sj + §9 ’ 


( 9 ) 


where k is the bulk modulus. In writing equation (9), we 
have neglected the extra terms which arise if there is a 
static elastic stress associated with the equilibrium state. 
It is easy to show that these terms are of order € yie id < 1 
relative to those retained. 

The crust is effectively a perfect electrical conductor 


at seismic frequencies, so the perturbed electric field is 
given by 

se= -~M xb - (io) 

Equation (10) combined with Maxwell’s equations, 


„ f „ 1 dSB 

VxSE = — - 

(ii) 

„ 4s- .. 1 dSE 

VxSB = — Sj + — 
c c at 

(12) 


may be used to relate 6j to £. 

The transmission coefficient for the energy flux between 
the crust and the magnetosphere should not be signifi- 
cantly affected by the polarization and propagation angle 
of the incident shear wave. We consider then the simplest 
possible case, a vertically propagating shear wave polar- 
ized such that £ oc z x B. For this special case V • £, 
6p y and 6p all vanish. Next, we take a harmonic time de- 
pendence, exp(— tu;*), for the perturbation variables and 
combine equations (7) - (12) to derive the wave equation, 

£( a s) + * w> * = 0 ' (13) 

where p and p are the effective shear modulus and density: 



Here cos a = B z jB. Equations (2) and (3) are used to 
obtain the dependence of p and p on depth. The magneto- 
elastic wave speed, v,, is equal to (p/p) 1 / 2 . In the magne- 
tosphere, z > 0, p = (B cos a) 2 /(47r), p — B 2 /(47rc 2 ), and 
the wave equation reduces to that for relativistic Alfven 
waves, 

j4 + = 0. (16) 

dz 2 c 2 

The transmission coefficient, T, is defined as the ratio 
of the transmitted to the incident energy flux. For a ho- 
mogeneous crust, it is easy to show that the transmission 
coefficient for a vertically propagating shear wave is given 
by the familiar formula 


4ZqZm 
(Zg + Zm) 2 ’ 

where Zc and Zm are the impedances of the crust and 
magnetosphere. Each impedance is the product of the rel- 
evant propagation speed and effective density. Thus, 

Z c = (p£) 1/2 and Z M = — . C ° g ~ • ( i8 ) 

47TC 

To determine the reflection and transmission amplitudes 
for an inhomogeneous crust, we perform two numerical in- 
tegrations of equation (13) from z <C 0 to z > 0 starting 
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from different lower boundary conditions. Then we choose 
the appropriate linear combination of the two boundary 
conditions such that the solution satisfies the radiation 
condition, = iwseca^/c, for z > 0. The curve in figure 
1 displays T as a function of frequency for various vertical 
magnetic field strengths. Only a small fraction T of the 
incident wave flux leaks into the magnetosphere. 



FREQUENCY (Hz) 


Figure 1: The transmission coefficient as a function of fre- 
quency for a vertically propagating wave. The upper and 
lower pairs of curves are for 10 12 and 10 11 G fields respec- 
tively. 


The fate of the reflected energy is clearly of interest. 
The angular momentum barrier (or equivalently, the 1 jr 
dependence of the horizontal wave number in a spherical 
star) and an increase in wave speed with depth both refract 
waves upward. Shear waves will not propagate in the fluid 
interior below the crust. Those which reach this region 
may be reflected if the rigidity drops sufficiently abruptly 
to zero at the inner boundary of the crust. Otherwise, as 
the waves slow down their radial wave vectors and am- 
plitudes will increase, and their energy will ultimately be 
dissipated as heat. 

The waves which do return will bounce many times off 
the surface and will spread throughout the crust. The 
characteristic storage time for wave energy in the crust, 
r(z/), is twice the time it takes the waves to cross the crust, 
2 t C y divided by T(u). For vertically propagating waves, 


other electrons, and free neutrons. We neglect chemical 
impurities because of the large uncertainty in their con- 
centration and type. At neutron drip, the damping time 
is of order 4 x 10 3 s for 10 4 Hz waves, the dominant scat- 
tering mechanism being due to phonons. This is probably 
an underestimate because of the likely presence of other 
sources of electron scattering, especially lattice imperfec- 
tions. Even so, it is still much longer than the typical 
timescale for a gamma-ray burst. 

Up to this point we have considered how a vertically 
propagating shear wave polarized orthogonal to the ambi- 
ent magnetic field couples to an Alfven wave. Since the 
magnetosphere can also support fast magnetosonic waves, 
it is natural to inquire whether more general seismic dis- 
turbances couple to them. Now, the dispersion relation for 
fast magnetosonic waves is isotropic and reads 

oj 2 — ( cfc ) 2 , ( 20 ) 

which differs from the anisotropic dispersion relation for 
Alfven waves, 


u> 2 = (ck cos x/>) 2 . (21) 

Here, k is the magnitude of the wave vector, and %j> is the 
angle between k and B. For fast magnetosonic waves and 
Alfven waves, Snell’s law reads 


sin B t c 

sin B{ v s 9 


( 22 ) 


sin Oj _ ccosV> 
sin B{ v a 

Here Bi and B t denote the angles of incidence and transmis- 
sion. At the reflecting layer, v M < c, and thus only waves 
having very small angles of incidence can couple to fast 
magnetosonic waves. For Alfven waves the only restriction 
is that the transmitted wave vector be almost orthogonal 
to B . 

The displacement amplitude of a propagating shear 
wave increases with height in the crust as a consequence of 
the conservation of energy flux. Because the strain ampli- 
tude, e, is proportional to the gradient of the displacement, 
and the magnitude of the wave vector, k = also in- 

creases with height, e reaches a maximum near the surface. 
The following argument relates the maximum strain am- 
plitude to the Alfven wave luminosity. 

We assume that the entire neutron star surface radi- 
ates Alfven waves. Then, the fractional perturbation of 
the surface magnetic field associated with luminosity L is 


t{u) ~ 


2 f ZQ dz 5 x 10" 6 |zg| 1 / 2 

T{u) J 0 T[y) 


(19) 


SB 

B 


— - ~ 6 x 10" 


/ £ X 1 / 2 ( B \ -1 

\10 38 ergs s _1 / \10 n G/ ^ ^ 


With the usual scalings, r ~ Is at = 10 4 Hz. 

Our analysis of the propagation of shear waves has 
neglected damping. Theoretical calculations of the degen- 
erate electron viscosity in a solid neutron star crust (Ref. 
9) included electron scattering by phonons, impurity ions, 


In terms of the surface displacement amplitude, equations 
(10) and (11) imply 


SB u> 
~B ~ c 


( 25 ) 
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Thus, wc estimate the total surface strain to be: 


6 ^^(l0 8 cms“ 1 ) (l0 38 ergss _1 ) (l0 11 G 

(26) 

With the chosen scalings, the maximum strain is danger- 
ously close to unity. A dynamic yield strain as large as 

0.1 is not unusual. Thus, we cannot predict whether shear 
waves transfer their energy to Alfven waves and/or crum- 
ble the crust, generating heat. However, it is clear that 
the neutron starquake model is not viable if bursts are sig- 
nificantly more distant, and consequently more luminous, 
than estimated above. 

4. ALFVEN WAVES IN THE MAGNETOSPHERE 

We have already discussed the magnetic field pertur- 
bations required for Alfven waves to carry the luminosity 
of a typical gamma-ray burst. Because Alfven waves trans- 
port energy without loss along the equilibrium field lines, 
their relative amplitudes vary according to 

(27) 

Equations (24) and (27) indicate that the Alfven waves are 
likely to become nonlinear not very far from the star. 

Up to this point, we have been making the implicit as- 
sumption that the plasma density in the magnetosphere is 
high enough so that the MHD limit applies to Alfven wave 
propagation. We shall now show that this is a questionable 
assumption. 

The corotation charge number density (Ref. 10), 



^^ 7 x 10 ’(io^g) (£) 1 


(28) 


sets a lower limit on the magnet ospheric plasma density, 
n. Here, P is the rotational period of the neutron star. 
Much higher densities may exist if accretion is occurring, 
or if an energetic process expells plasma from the stellar 
atmosphere. 

For low plasma densities a nonlinearity may occur 
even at relative field amplitudes SB / B < 1. Because the 
wave vector is nearly orthogonal to the magnetic field near 
the surface, equation (10) implies that the wave vector is 
nearly parallel to SE. Therefore, the displacement current 
and V x SB are nearly orthogonal. From equation (12), 
we see that a substantial physical current is necessary to 
support the wave. If the plasma density is so low that 
the required drift velocities exceed c, the wave is charge 
starved. This occurs where 


SB 2 nec 

~B ~ i/o B 



Based on the values of SB/B given by equation (24), the 
Alfven waves are expected to be charge starved close to 
the stellar surface unless n > n cr . The fate of such waves 
requires further investigation. 

Once the Alfven waves go nonlinear, a substantial part 


of the wave energy is probably transferred to the electrons. 
An upper limit to the Lorentz factors is obtained by bal- 
ancing the electrostatic acceleration with curvature radia- 
tion reaction. We find 



5. CONCLUSIONS 

Starquakes have the following virtues. They can easily 
provide the energy necessary to power a single burst. They 
release energy in the low entropy form of seismic waves 
which are transformed into Alfven waves on timescales 
characteristic of observed gamma-ray bursts. The Alfven 
waves may become nonlinear far from the the stellar sur- 
face and accelerate the ambient plasma to energies that are 
limited by radiation reaction losses due to the emission of 
gamma-rays. Thus, it is possible to imagine an efficient 
conversion of elastic energy to gamma-ray energy. 

Of course, the above scenario is replete with uncer- 
tainty. It is quite possible that plastic flow, and not brittle 
fracture, is the manner in which stress is relieved in neu- 
tron star crusts. The source of the free energy necessary to 
replenish the crustal stress in old neutron stars is not iden- 
tified. The conversion of seismic waves into Alfven waves, 
on the appropriate timescales and without the production 
of excessive crustal strains, requires that magnetic fields 
> 10 n G cover at least part of the stellar surface. The 
propagation of Alfven waves in low density, neutron star 
magnetospheres is poorly understood. Their behavior and 
the manner in which they accelerate particles probably de- 
pends upon the degree to which the neutral plasma density 
exceeds the corotation change density. 
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